Effects of optical fields on the tunneling time of chiral electrons in graphene 
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The influences of optical fields on the tunneling time in graphene are investigated in real time 
using the finite-difference time-domain method. The tunneling time of electrons irradiated by an 
optical field is significantly different from that observed in traditional quantum tunneling. We 
found that when the barrier width increases, the group delay becomes constant for the reflected 
wave packet, but increases linearly for the transmitted wave packet. This peculiar tunneling effect 
can be attributed to current leakage in a time-dependent barrier generated via the optical Stark 
effect. 

PACS numbers: 42.65.-k, 68.65.-k, 73.40.Gk 
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Quantum tunneling time has received much attention 
since MacColl pointed out that it takes no approximate 
time [1.-3J . The question about how long it takes an elec- 
tron to tunnel through a potential barrier is still replete 
with controversy. The debates center around the defini- 
tion of tunneling time and its exact physical meaning. 
Hartman calculated the group delay or phase time and 
found that the the group delay Tg becomes constant while 
the barrier length increases U]. Thus, with a wider bar- 
rier, superluminal group velocities can be observed. Re- 
cently, Winful proposed that the group delay in tunneling 
represents a lifetime of stored energy escaping through 
both sides of the barrier, not a transit time |3|, |5|, |6|- 
However, unlimited group velocity is not a meaningful 
concept in tunneling and does not imply superluminal- 
ity. 

In the meantime, many optical and acoustic experi- 
ments have been carried out to determine tunneling time. 
Steinberg et al. [7| measured the time delay for a pho- 
ton to tunnel across a one-dimensional photonic band- 
gap material. They reported that the group velocity for 
single-photon tunneling is about 1.7co, where cq is the 
speed of light in vacuum. The authors attributed super- 
luminality to the fact that the wavepacket may be re- 
shaped in the tunneling process. Longhi et al. |8| used a 
relatively long optical pulse (380ps) and conformed that 
there is no distortion in the the tunneling process. Even 
though the exact physical meaning of these experimental 
results is controversial, there is no doubt that there is a fi- 
nite duration for the photon or phonon tunneling process. 
However, few direct experiments have shown that there 
is finite tunneling time in the quantum particle tunneling 
process. 

There are some difficulties in the measurement of the 
tunneling time of electrons. For instance, the coherence 
time of electrons must be long enough to ensure that the 
tunneling process is an elastic transport. It seems that 



the graphene is an ideal candidate material for this. The 
reported electron mobility in graphene at room temper- 
ature is in excess of 15,000 cm^V~^s~^, and the limit of 
electron mobility at room temperature is about 200,000 
cm'^V~^s~^ in theory [3-'ll'|. To measure the tunneling 
time, a high time-resolved technique must be used, so the 
injected electrons must be generated using ultrafast laser 
beams (e.g., generated via coherent one- and two-photon 
or the barrier must be controlled 
Furthermore, be- 



15, 1 



absorption [l2| - [l3 |) 
by the ultrafast laser beams 
cause of the existence of zero-point fields, the influences 
of electromagnetic fields on the tunneling time must be 
carefully treated. When the electromagnetic fields are 
included, the barrier is time dependent. The situation is 
quite different when a time-dependent potential barrier is 
taken into account. For instance, for an opaque rectangu- 
lar barrier with a small time-dependent modulation, the 
traversal time of Biittiker and Landauer is proportional 
to the barrier width D |17| . 

In this Letter, we investigate the influences of electro- 
magnetic fields on the tunneling time of Dirac electrons 
in graphene. We consider a rectangular potential bar- 
rier with height Vq and width D [see Fig. 1 (a)]. The 
Fermi level lies in the valence band in the barrier and in 
the conduction band outside the barrier. A F-direction 
polarized laser beam is propagated perpendicular to the 
layer surface with a detuning Aq — Ei, ~ huj. We choose 
Aq > to ensure that there is no interband absorption 
inside the barrier, and Huj <^ 2Ek to guarantee that the 
influence of the optical fleld outside the barrier can be 
neglected. In a single layer graphene, a perfect transmis- 
sion through a potential barrier in the normal direction is 
expected, and there is no Hartman effect |18l421| . How- 
ever, if the barrier is irradiated by a intense optical fleld, 
the conduction band and valence band is mixed, the chi- 
ral symmetry of the Dirac electrons is broken, and the 
perfect tunneling is strongly suppressed 16|. Thus, the 
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FIG. 1: (color online), (a) Schematic of the scattering of 
Dirac electrons by a square potential irradiated by an elec- 
tromagnetic field. Ba, Bin, and Bout denote the absorbing 
boundary, incident boundary, and output boundary, respec- 
tively, (b) Schematic of the one-dimensional Yee lattice in 
graphene. (c) The evolution of a wave packet through a bar- 
rier with pump intensity I^, = 20 MW/cm^, Ao — 5meV, 
D — 150 nm, and Vo ~ 400 meV. The light grey shows the 
barrier area. 



Hartman effect may appear if the electromagnetic field is 
included. 

In order to study such a time-dependent scattering 
process, we employ the finite-difference time-domain 
(FDTD) method to solve the time-dependent Dirac equa- 
tion numerically. In the FDTD method, the time- 
dependent Dirac equation in K point is replaced by a 
finite set of finite differential equations 
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where (i, k) — (iAx, kAt) denotes the grid of point of 
the space [see Fig .1(b)], *(r,t) = [CA(r, t), CB(r, t)] is 
the wave function of Dirac electrons, vp ~ lO^m/s is 
the Fermi velocity, Vb(r) is the height of the potential 
barrier, and V12 and V21 are the matrix elements of the 
interaction Hamiltonian Hint 



m 



B.i.nt = -hevF{A^(J^+Ayay) = h[ y q^ ) > (2) 



0.12 



I 0.06 
Q 

0.03 



(a) ... 
. / 

P 



7 J -.-l=20MW/crrf 
-"-l=30MW/cm^ 
l=50MW/cm' 
100 200 300 400 

D (nm) 



1.5 
1.0 
0.5 
0.0 



/, Y.-l=20MW/cm 
yPl - |=30MW/cm' 
-'-l=50MW/cmS«' 



300 600 900 1200 

D (nm) 



0.0020 r 




00 400 600 800 10001200 

X(nm) 



FIG. 2: (color online), (a) Group delay for the reflected wave 
packet and (b) group delay for the transmitted wave packet 
as a function of barrier width D at different pump intensities 
for Ao = 5 meV. (c) Snapshots of the spatial probability den- 
sity distributions on the left side of the barrier for different 
pump intensities. The solid lines are the best fit using ex- 
ponential functions, (d) Snapshots of the spatial probability 
density distribution on the right side of the barrier for I^^ = 50 
MW/cm^ 



where e is the electron charge and {A^TAy) are the vector 
potentials of the electromagnetic field. 

Thus, by numerically solving Eq. P^ and Eq. (ITb)) di- 
rectly in the time domain, we can demonstrate the prop- 
agation of a wave packet through a barrier in real time. 
For computational stability, we choose the space incre- 
ment Ax = 0.1 nm and the time increment At = 5xl0~^ 
ps. As an example, the evolution of a wave packet 
through a barrier is shown in Fig. 1(c). At the in- 
put boundary i3i„, a Gaussian electronic wave packet 



Ca^ Cb ^ exp 



47r(t-to)^ 



is injected. For the conve- 



nience of demonstration, a short pulse is used: the peak 
position to — 1.2 ps, and the pulse width r — 0.8 ps. 
From Fig. 1(c), we find that when the sample is irra- 
diated by an intense nonresonant laser beam, a reflected 
wave packet appears, and the perfect transmission is sup- 
pressed. 

Under an intense optical field, the light-induced band 
shift can also create a dynamic gap at small detuning. 
As is seen in traditional quantum tunneling (e.g., the 
Hartman effect), the tunneling time of Dirac electrons 
should be constant while the barrier width increases. To 
verify this, we studied the group delay (i.e., the delays 
of the peak of the reflected and transmitted pulse) of 
Dirac electrons through a potential barrier generated by 
an intense light beam. In order to reduce the distortion, 
a relatively long plus is used: the peak position tg = 5.0 
ps, and the pulse width r = 3.3 ps. 

Fig. 2(a) shows the group delay for the reflected wave 
packet as a function of barrier width D for different pump 



intensities. As seen in traditional quantum tunneling, the 
group delay is saturated by increasing the barrier length. 
For I^ = 20 MW/cm2, I^ = 30 MW/crn^, and I^ ^ 50 
MW/cm^, the corresponding saturated delays are about 
0.112, 0.0893, and 0.0662 ps, respectively. The delay 
can also be explained by the tunneling depth or dwell 
time. As shown in Fig. 2(c), the spatial probability den- 
sity distributions in the barrier are well fitted by using 
the exponential function F — A^e^^/^'^. For I^^ — 20 
MW/cm^, /„ = 30 MW/cm2, and /^ = 50 MW/cm^, 
the corresponding tunneling depths xq are 55, 44, and 33 
nm, respectively. In this way, we extract that the corre- 
sponding tunneling delays io = '2'Xo/vf are about 0.11, 
0.088, and 0.066 ps, respectively. These results are con- 
sistent with the saturated delay obtained in Fig. 2(a). 
Another interesting phenomenon is that the group delay 
for the reflected wave packet takes on the quantized val- 
ues T = nTi {n = 1,2,3...). We find that Ti « 0.011 ps 
is equal to the period of the optical field. The quantized 
group delay is therefore caused by the optical modulation 
of the reflected wave packet. 

The case is quite different for the transmitted wave 
packet. As shown in Fig. 2(b), when the barrier width 
D > 900 nm, the group delay increases linearly with in- 
creasing barrier width. The group velocity is about 10^ 
m/s and is the same as the Fermi velocity of Dirac elec- 
trons in graphene. This result can be explained by the 
time-dependent potential barrier. In a time-dependent 
potential barrier, the energy storage depends on time. 
Variations in energy storage will lead to an extra leakage 
current. The dynamic gap caused by the optical field is 
therefore not a complete gap, and the probability density 
distribution on the right side of the barrier is no longer 
exponential decay [see Fig. 2(d)]. Since the amplitude of 
the extra leakage current is quite small, it has little ef- 
fect on the reflected wave packet but still determines the 
group delay for the transmitted wave packet at large bar- 
rier widths. For small barrier widths, since the tunneling 
current is much larger than the extra leakage current, 
the time delays of the reflected and transmitted wave 
packets are equal. From Fig. 2(b), we see that the am- 
plitude of the extra leakage current strongly depends on 
the pump intensity. Specificially, for I^ = 50 MW/cm^, 
the group delay increases linearly with increasing barrier 
width when D > 700 nm, but for /„ = 20 MW/cm^ 
the linearly increasing delay appears when D > 900 nm. 
However, the group velocity of the extra leakage current 
is independent of pump intensity. Thus, even with a quite 
small time-dependent modulation, no "Hartman effect" 
occurs. 

The existence of current leakage in a time-dependent 
barrier can also be confirmed by the shape of tunneled 
pulses. In a thin barrier, the amplitude of the tunneling 
current is much larger than the amplitude of the extra 
leakage current induced by the time-dependent modula- 
tion. The tunneling time of the transmitted and refiected 
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FIG. 3: (color online). Incident (black lines), tunneled (green 
lines), and reflected (red lines) pulses with pump intensity 
/„ = 30 MW/cm^ for the diflterent barrier width (a) D = 100 
nm, (b) D = 500 nm, and (c) D = 900 nm. The insets 
shows the normalized tunneled pulse overlaid with the inci- 
dent pulse. 



wave packet is equal [see Fig. 3(a)], and the distortion is 
quite small [see the inset of Fig. 3(a)]. When the thick- 
ness of the barrier increases, the tunneling rates decrease 
rapidly. In a thick barrier, the amplitude of the extra 
leakage current is comparable with that of the tunneling 
current. Since the tunneling times of the extra leakage 
current and the tunneling current are different, a seri- 
ous distortion can be found in the tunneled pulse [see 
Fig. 3(b)]. However, if the width of the barrier is large 
enough (e.g., D = 900 nm), the tunneling rate is very 
small, and the extra leakage current is the main contrib- 
utor to the tunneled pulse. Thus, the tunneling times of 
the transmitted and reflected wave packets are different 
[see Fig. 3(c)], and there is no distortion on this scale 
[see the inset of Fig. 3(c)]. 

The tunneling delay for a time-dependent potential 
barrier might help us to understand quantum tunneling. 
For example, traditional definitions of tunneling time are 
unsuitable for the system we studied. In the traditional 
definitions, the tunneling time of the transmitted and 
refiected wave packets are equal for a symmetrical bar- 
rier. More importantly, if the quantum fiuctuation or the 
zero-point field is included, all potential barriers are time 
dependent. 

In conclusion, we have calculated the influence of op- 
tical fields of chiral tunneling time in graphene using the 
FDTD method. We find that the group delay of the 
refiected packet is also saturated as the barrier width in- 
creases. However, the delay increases linearly with bar- 
rier length for the transmitted wave packet. This peculiar 
tunneling effect is attributed to current leakage in a time- 



dependent barrier generated via the optical field. Thus, 
the zero-point field may have an important influence on 
the tunneling time of electrons, and should be carefully 
treated. These unique transport properties of dressed 
Dirac electrons in graphene might be important in the 
understanding of quantum tunneling. 
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